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ABSTRACT 

The Gamma-Ray Burst (GRB) prompt emission is believed to be from 
highly relativistic electrons accelerated in relativistic shocks. From the GRB 
high-energy power-law spectral indices (3 observed by the Burst and Transient 
Source Experiment (BATSE) Large Area Detectors (LAD), we determine the 
spectral index, p, of the electrons’ energy distribution. Both the theoretical 
calculations and numerical simulations of the particle acceleration in relativis¬ 
tic shocks show that p has a universal value « 2.2 — 2.3. We show that the 
observed distribution of p during GRBs is not consistent with a (^-function 
distribution or an universal p value, with the width of the distribution > 
0.54. The distributions of p during X-ray afterglows are also investigated and 
found to be inconsistent with a (5-function distribution. The p-distributions in 
blazars and pulsar wind nebulae are also broad, inconsistent with a (5-function 
distribution. 
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1 INTRODUCTION 


GRBs are observed to have non-thermal spectra during its prompt emission phase (Band et 
ah 1993). It is widely believed that the synchrotron radiation and/or the inverse Compton 
scattering are the likely emission mechanism(s) for CRB’s prompt hard X-ray and 7 ray 
emission. The electrons accounting for these emissions are thought to be accelerated in 
relativistic shocks in CRBs. According to the shock diffusive acceleration model, particles 
are accelerated when they repeatedly cross a shock front, and the competition between the 
particle’s energy again and escape probability per shock crossing cycle leads to a power-law 
spectrum for the particles: 

A^( 7 )(i 7 oc 7 “^(i 7 , (1) 

where 7 is the Lorentz factor of the particle (e.g., Blandford & Ostriker 1978). For non- 
relativistic shocks, the value of p depends on the the compression ratio of the flow stream 
across the shock; while in relativistic or ultra-relativistic shocks, which is most likely the case 
in CRBs, analytical and numerical studies show that p has an “universal” value, ~ 2.2 — 2.3 
(Kirk et al. 2000, Achterberg et ah 2001, Bednarz & Ostrowski 1998, Lemonine & Pelletier 
2003). 

The work reported in this paper is to investigate the “universality” of the power-law 
index p for CRBs, which we calculate directly from the high-energy (0.1 — 2 MeV) photon 
spectrum of CRBs (Preece et al. 1998, 2000), assuming the spectrum is from synchrotron or 
synchrotron self-inverse-Compton emission of the power-law distributed highly relativistic 
electrons, using the relations between p and the spectral index , (3, of the high-energy power- 
law photon spectrum. 

In §2, we describe the CRB spectral data set used and the process of determining the 
parent p-distribution. In §3, we examine the contribntions from the spectral £t procednre and 
the time averaging effect to the dispersion of the parent distribntion of p. The p-distribntions 
derived from BeppoSAX CRBs and from HETE-2 CRBs, X-ray flashes and X-ray rich CRBs 
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Electron spectral indices in GRBs & other relativistic sources 3 


are presented in §4. We determine the p-distributions for X-ray afterglows in §5 and for 
blazars and pulsar wind nebulae in §6. A summary and discussions are given in §7. 


2 THE DISTRIBUTION OF P IN GRBS 
2.1 The GRB spectral sample 


For our analysis, we use the BASTE GRB Spectral Catalog presented by Preece et ah 
(2000). In the catalog, the time sequences of spectral fit parameters for 156 bright bursts 
are presented, using mostly the high energy and time resolution data from the Large Area 
Detectors (LAD), which covers an energy range of typically 28 - 1800 keV. All bursts have 
at least 8 spectra in excess of 45 a above background. The spectral models used in fit are (i) 
‘Band’ function; (ii) Comptonized spectral model (a power-law with an exponential cut-off); 
(hi) Broken Power-Law model; and (iv) Smoothly Broken Power-Law model. The ‘Band’ 
function, the one used most frequently, is an empirical function (Band et ah 1993) 


N{E) = A 


(E/100)"exp[-E(2 -h a)/Epeak], 


^ ^ 2+o-^PeaA: 


L 100(2+o) 


exp(/3 - a)(E/100)^, E > ^Ep^ak 
where N{E) is the photon counts, A is the amplitude, a is the low-energy spectral index, 
f3 the high-energy spectral index, and Epeak is the peak energy in the vEy spectrum (when 

/9<-2). 


Since we are here caring about the high-energy power-law portion of the GRB spectra, 
and also because one possible source of systematic error in the spectral parameter determi¬ 
nation arises in selecting different spectral models for different bursts (Preece et al. 2002), 
only the spectral parameters of those ’Band’ function fitted spectra are selected for our 
analysis. 

One of our major concerns is to select the sample of spectra for which (3 is reliably 
determined. The BATSE burst signal-to-noise ratio decreases at higher energies as a result 
of fewer photon flux and the decreased detector efficiency. In particular, (3 may not be well 
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determined if -Epeafc is close to the higher limit of the LAD energy range, E^ax (~ 2 MeV) 
(Preece et ah 1998), thus we must choose those spectra with Epeak much lower than Emax- 
Therefore, we select the spectrum for which 100 keV < Epeak < 200 keV and the error in (3 
is less than 0.1 |/5|. This gives a total sample of 395 spectra for 78 bursts. 


2.2 Distribution of p and its narrowing 

For electrons’ distribution given by a power-law: 

^(7e) oc 7"^, for 

'Je ^ 'ymint ( 2 ) 

the emergent high energy synchrotron spectrum is asymptotically a power law function: 
F7 oc for Um < n < Ue (“slow cooling” regime) and oc for v > Ve (“fast cool¬ 

ing” regime), where = t^synilmin) is the synchrotron injection frequency, and = nsyn{'jc) 
is the synchrotron cooling frequency above which the synchrotron energy loss becomes im¬ 
portant. 

The spectral index, p, of shock accelerated electrons is associated with the high-energy 
power-law photon index, (3, of GRB photon spectrum, by either f3 = —p/2 — 1 (“fast cooling” 
regime) oi j3 = —{p + l)/2 (“slow cooling” regime) depending on relative positions of Vm and 
Vc and on which portion of the spectrum is detected. There is one regime, z/^ < z/ < z/^, m 
which f3 = —3/2, independent on p. This case can be ruled out by discarding those spectra 
with (3 > —3/2 from our sample. We found only one with (3 > —3/2 in the BATSE sample 
of 395 spectra and discarded it. 

Piran (2004) argues that the fast cooling must take place during the GRB prompt phase 
and the reasons are: (i) the relativistic shocks must radiate their energy efficiently, to avoid 
a serious inefficiency problem; (ii) the electrons must cool rapidly in order that the fast 
variability could be observed. But there is no hrm evidence to date that could rule out the 
slow cooling case for the GRB itself, since it is difficult to measure the values of 7 c and 
7 mm for a specihc burst. Thus in our analysis, we assume that each GRB spectrum above 
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Electron spectral indices in GRBs & other relativistic sources 5 
Epeak could bc lu either slow cooling or fast cooling regime, so as to minimize the width of 
p distribution. 

First we plot distribution of p by assuming all spectra are in fast cooling regime. Then we 
make the distribution narrower by relaxing this constraint. Basically the narrowing process 
is to move some left-hand part of the distribution to the right by adding 1 to p and assuming 
this part of sample are in slow cooling regime, since there is a difference of 1 about p value 
between the two regimes. The algorithm used is described below. 

Several algorithms are implemented to get the narrowest distribution. In the most straight¬ 
forward one, each spectra has the freedom of calculating p from (3 either in “fast cooling” 
or “slow cooling” regime, so the number of possible distributions is 2^ for a sample of N 
spectra. The distribution having the smallest standard deviation is chosen as the narrowest 
one. This algorithm works well only for iV < 20 because of the computer program’ run¬ 
ning time. For N > 20, we divide the overall range of the sample’s (3 distribution into 20 
equal-width bins and treat the spectra with /5’s located in each bin indistinguishably. Then 
we apply the first algorithm to the 20 bins. In an alternative algorithm, we start with the 
histogram of p calculated in the “fast cooling” regime and mark a demarcation line within 
and close to the lower limit of the range of p. Then all p’s at left to the line in the his¬ 
togram are moved to right by adding 1 to p, and the new histogram’s standard deviation is 
calculated. Repeat it after shifting the demarcation line rightward by a step of 0.01 on the 
p-axis. Finally the smallest standard deviation, hence the narrowest distribution, is found. 
It turns out that both algorithms give the same results for most of the samples presented 
in this paper. For one sample where minor difference exists between two algorithms’ results, 
we use the narrower one. 

We show the results of the analysis for BATSE bursts in Figure 1. Note that all the errors 
presented in this paper are at 1 a level. The parent distribution of p for BATSE bursts has 
a width of 0.54 at a 14-a confidence level. The method that estimates the mean and the 
width of the parent distribution of p is described below. Note that the mean value of p is 
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~ 3, substantially larger than that for the distribution before the minimization, which is an 
artifact of choosing some of the spectra to be in the “slow cooling” regime, equivalent to 
moving the left part of the histograms in the upper panels rightward, in order to minimize 
the width of the distribution. 


2.3 Statistical description of the narrowness of p’s distribution 

The observed distribution of p plotted in Fig. 1 is a convolution of the measurement error 
distribution and the true distribution (or parent distribution) of p’s. What we want to know 
is the true distribution p. We use the maximum likelihood method to estimate the true 
p-distribution. Let us say the true distribution of p is Gaussian, 


P(p) = 


27T(!T^ 


■expl-lk-^’L 




(3) 


Further, we assume the measurement errors have Gaussian distributions too. Then the prob¬ 
ability distribution for any one measurement (p*, Ui) is the convolution of two Gaussians, 
which is the Gaussian 


P(p*,ai,p,ap) = 

The likelihood function for the set of n measurements pi, ai is 


1 {Pi - Pf 
'2 + 


L = I[ 




exp[— 


1 {pi - pf 

2 


(4) 


(5) 


The principle of the Maximum Likelihood Estimate is that, the best estimates of p and 
are the ones that maximize L. Take 


1 ^ (pi - p)2 1 ^ 

I = lnL = 2 - T ^ 

^ i ^ i 

then the maximum occurs when the following equations 
dl 


dp 


= 0 , 


dl 




= 0 


( 6 ) 

(7) 

( 8 ) 
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have their solution at p = p and Up = Up, where symbolizes the best estimation of the 
parameters. If we assume that the distribution of p and Up are both Gaussian, then one can 
show that the variances of p and Up are 


2 

■ OH 

■ 

^p = - 

gp2 



1-1 


a. 2 = 

Up 


OR 





1 -1 


(9) 

( 10 ) 


respectively. So the best estimate of the parameters of true distribution of p are obtained by 
nnmerically solving eqnations (7) and (8), and their associated errors are calculated through 
equations (9) and (10). 


3 SYSTEMATIC ERRORS IN (3 

3.1 The ‘Band’ function fit to the spectra 

Preece et ah (2000) carried out a Band fnnction £t to GRB spectra observed by BATSE, 
and this way determine the high energy power-law index (/?) and the random error in (3 due 
to error in the observed spectral energy distribntion. There is also a systematic error in (3 
resulting from the finite bandwidth of the BATSE detector, which was not reported in Preece 
et ah, and we estimate it here. The purpose of this exercise is to estimate the contribution 
of this systematic error, and its dependence on the peak of the spectrum (Epeak), to the 
dispersion in the p-distribution. 

The systematic error arises because the synchrotron spectrnm does not make a sharp 
transition from one power-law index to another when one crosses a characteristic freqnency. 
In particular, the steepening of the spectrum to above the synchrotron and cooling 

frequencies does not occur suddenly at Epeak, but instead the spectrum approaches this 
theoretical value asymptotically at E 3> Epeak- 

Since the spectrum is observed in a finite energy range, the measnred spectral index will 
always be somewhat smaller than the trne asymptotic value by an amonnt that depends on 
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the ratio of E^ax and Ep^ak {Emax is the highest energy photon that the detector is sensitive 
to). The larger the Emax/Epeak is the smaller the systematic error in f3 would be, and this 
dependence on Epeak causes some broadening of the observed f3 distribution. 

To estimate this systematic error we generate synthetic spectra with different values for 
Epeak^ and carry out a Band function £t to the synthetic spectra to determine f3 and its 
deviation from the true asymptotic value. 

The synthetic synchrotron spectra is calculated for a relativistic homogeneous shell. The 
electron distribution function behind the shock is taken to be a single power-law function: 
iV( 7 e) oc 7 “^, for 7 > 7min, where meC^'Jmin is the minimum electron energy after they cross 
the shock front. The magnetic held in the shell is taken to be uniform and the energy density 
is the held is some fraction (e^) of the thermal energy density of the shocked fluid; and 
eB are chosen so that the peak of the spectrum, Epeak, is at some desired value. As electrons 
move down-stream from the shock front they cool via the synchrotron and inverse-Compton 
processes, and their distribution function is modified. We calculate the effect of this cooling 
on electron distribution functions using a self consistent scheme described in Panaitescu & 
Meszaros (2000) and McMahon et al. (2006). 

The synchrotron spectrum, for a given electron distribution, in the shell comoving frame 
is calculated as described in detail by Sari et al. (1998) (also see section 2.2). The spectrum 
in the observer frame is calculated by integrating the spectral emissivity in the comoving 
frame over the equal-arrival-time surface as described in Kumar & Panaitescu (2000). Errors 
are then added to this spectrum in a way that mimics the real GRB spectrum. 

The synthetic spectrum for a known p is fitted to the Band function in a finite energy 
range corresponding to the BATSE energy coverage. By varying Epeak of the generated 
spectra we determine the discrepancy between fitted value and “true” value of /? as a function 
of Emax/Epeak- The results are shown in Figure 2. We find the fit always gives a smaller f3 
(in absolute value) than the true asymptotic value and that the “observed” f3 does indeed 
depend on Emax- The error in f3 is about 10% when Emax/Epeak is order unity, whereas 
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the error is 5% when E^ax/Epeak 20. The error also depends on the p value as shown 
in Figure 2; for Epeak located between 100 keV to 200 keV, Emax= 1-8 MeV, and p= 2.5, 
the contribution of this systematic error to the dispersion in (3 is less than 1.3% - the 
corresponding contribution to the dispersion in p is dp < 0.03. 

We have also carried out a similar calculation for the synchrotron self- inverse-Compton 
(SSC) spectrum for a population of synchrotron electrons. The incident photons are the 
synchrotron photons due to the same population of electrons that contribute to inverse- 
Compton scatterings. The synchrotron radiation is taken to be homogeneous and isotropic 
in the shell comoving frame, and its spectrum is calculated as described above. The overall 
SSC spectrum is obtained by the convolution of the synchrotron spectrum and electron 
energy distribution using equation (7.28 a) in Rybicki & Lightman (1979). The curvature 
in the SSC spectrum is due to the convolution of the incident spectrum and the electron 
distribution, and we hnd that the asymptotic value for the SSC power-law index is reached 
when E^ax!Epeak 100. For this reason we hnd that for the SSC case, the systematic error 
in is ~ 13% for the typical E^ax/Epeak in BATSE bursts. The dispersion in p caused by 
Epeak being distributed between Epeak= 100 keV to Epeak= 200 keV is, however, small - 
(Tp < 0.04. 

These results show that the discrepancy between the htted value and the “true” value 
of f3 is small and dependent on Emax/Epeak, but its dependence on Emax/Epeak is too small 
to account for the observed dispersion in the p distribution. 


3.2 Time-averaging effect 

Another source of systematic error in (3 is the time-averaging of multiple spectra undergoing 
spectral evolution, i. e., Epeak evolving with hux (Ford et al. 1995, Crider et al. 1999). The 
hux-weighted time-averaging of multiple ‘Band’ spectra may distort the intrinsic high-energy 
power law. 

To examine this effect, we select BATSE time-resolved spectra with Epeak in 100 - 200 
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keV and in 200 - 300 keV, respectively, divide them into non-evolving gronps and evolving 
gronps, and analyze their p distribntions separately. The results are shown in Table 1. We 
hnd the evolving spectra groups tend to have flatter p or (3, which may be an outcome the 
time-averaging effect. But the widths of p distributions for two groups are consistent with 
each other, showing that the time-averaging does not contribute to observed dispersion in p 
in Figure 1. 

The time-averaging effect is further examined when we use an early BATSE spectral 
catalog by Band et al. (1993) in which the time-integrated spectrum of each burst is fitted 
with the ‘Band’ function. We restrict our samples to those with Ep^ak < 300 keV, and error 
in (3 less than 0.1 \j3\, which gives a sample of 32 spectra from the catalog of 54 GRBs. 
The p distribution is shown in Figure 3. Comparing with Figure 1, one can see that it has 
approximately the same cjp as that for the time-resolved GRB spectra. This supports that 
the time-averaging effect has no impact on the observed dispersion in p. 

4 P-DISTRIBUTIONS FOR BEPPOSAX GRBS AND HETE-2 XRFS, 

XRRS AND GRBS 

We also analyzed a sample of 11 GRBs observed by BeppoSAX. The combined (2 - 700 
keV) Wide Field Cameras (WFC) and Gamma-Ray Burst Monitor (GRBM) spectra for 
these bursts are fitted with the ‘Band’ function by Amati et al. (2002). The narrowest 
distribution of p for this sample is shown in Figure 4 left panel. It has the same estimated 
mean value of p as in the BATSE bursts, and the width of the parent distribution for p is 
consistent with that for the BATSE bursts. The larger errors in < p > and ap are due to 
the smaller size of the BeppoSAX sample. 

Sakamoto et al. (2005) present a catalog of X-ray flashes (XRFs), X-ray-rich (XRR) 
GRBs and GRBs observed by HETE-2 WXM (2 - 25 keV) and FREGATE (7- 400 kev) 
instruments. Among 45 bursts in the catalog, 16 bursts have measured high-energy power- 
law photon index, (3, which is obtained through the spectral fit with the ‘Band’ function or 
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a single power-law model. For those XRF spectra fitted by a single power law, it is found 
that (3 <-2. Sakamoto et ah (2005) explain this as that we are observing the high-energy 
power-law portion of their “Band”-function spectra. Two GRBs (GRB 020813 and 030519) 
for which “Band” model is used have Ep^ak lying near or above the upper limit of FREGATE 
energy range, so we exclude them here. We also exclude XRF 030528 which has a large error 
in j3. The hnal HETE-2 sample we considered comprises 7 XRFs, 4 XRRs and 2 GRBs. The 
p distribution is shown in Figure 4 right panel. 


5 P-DISTRIBUTION FOR X-RAY AFTERGLOWS 

We also determine the distribution of p during the X-ray afterglows. We use a catalog of X- 
ray afterglows observed by BeppoSAX compiled by De Pasquale et ah (2005) and a catalog 
of X-ray afterglows observed by Swift (O’Brien et ah 2006). In De Pasquale et al. (2005) ’s 
catalog, 15 X-ray afterglow spectra are htted with a Galactic-and-extragalactic absorbed 
single power law. We use 14 out of them for our analysis and exclude GRB 000210 which 
has an extremely large error in measured j3. In O’Brien et al. (2006) ’s Swift catalog of 40 
X-ray afterglows, we select samples with small errors, a(/5j) < 0.1|/5i|, and discard a sample 
with extremely large |/3| (= 5.5). We also discard 4 samples with jdi > —3/2 because these 
(3 values indicate the X-ray band probably lies between and (^'c < < ^m)i where 

the asymptotic spectral index is j3 = —3/2 and carries no information about p. This gives 
28 samples from the catalog. 

The p-distributions for the two afterglow samples are shown in Figure 5. For the Bep¬ 
poSAX afterglows, the narrowest distribution is consistent with a h-function distribution 
within 1 a errors; for the Swift afterglows, it is not. The smaller estimated width of the 
parent p-distribution for BeppoSAX afterglows, we suspect, is due to larger errors in photon 
indices Pi of the BeppoSAX sample, < apP) >= 0.26, than the Swift sample which has 
< aPP) > = 0 . 10 . 
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6 DISTRIBUTION OF P IN BLAZARS AND PULSAR WIND NEBULAE 

6.1 Blazars 

Blazars are active galactic nuclei with the relativistic jet pointed toward us. The nonthermal 
spectra of blazars are due to synchrotron or/and inverse Compton emission of relativistic 
electrons accelerated by shocks within the jet (Blandford & Konigl 1979, Sikora et ah 1994). 

Donato et ah (2005) present a spectral catalog of BeppoSAX six years of observations of 
Blazars at 0.1 - 50 keV. This catalog comprises three classes of blazars, namely low-luminosity 
sources (High-energy peaked BL Lacs, or HBLs), mid-luminosity sources (Low-energy peaked 
BL Lacs, or LBLs) and high-luminosity sources (Flat Spectrum Radio Quasars, or FSRQs). 
The three classes have different locations of synchrotron peak. X-rays from HBLs are likely 
to be above the peak of synchrotron spectrum, thus have steep X-ray spectra {{3 < -2), while 
FSRQs and LBLs in X-ray band have more contribution from inverse Compton component 
and thus have flatter spectra. 

From this catalog we use 44 spectra of 33 HBLs (some sources have multi-epoch spectra) 
that are best htted by single power-laws. The errors of htted photon indices reported in 
Donato et al. (2005) are at 90% conhdence level which we convert to 1-a errors. The distri¬ 
bution of p derived from their photon spectral indices is shown in Figure 6. We hud that the 
distribution of p for blazars is not consistent with a (5-function distribution: Cp = 0.22 ±0.03 
after the narrowing. 

6.2 Pulsar wind nebulae 

Power-law nonthermal spectra are also often observed in pulsar wind nebulae (PWNs) of 
rotation-powered pulsars. The nebular emission is the synchrotron radiation from charged 
particles heated by the termination shock in relativistic outflow (winds) from the pulsar 
(see Arons (2002) for a review). Gotthelf (2003) presents a catalog of nine bright Crab-like 
pulsar systems with Chandra observations and the photon indices of pulsar nebulae, Ppwn, 
and their 90% conhdence errors are provided. We derive the distribution of p from Ppwn 
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with the Ppwn errors converted into 1 a errors and hnd that dp = 0.59 ± 0.15, < p >= 
1.72 ± 0.20 assuming the X-ray band is in the fast cooling regime. After narrowing, the 
narrowest distribution has dp = 0.24 ± 0.07, < p >= 2.04 ± 0.09. 

7 SUMMARY AND DISCUSSIONS 

Motivated by theoretical calculations and numerical simulations showing that the shock- 
accelerated electrons in relativistic shocks have a power-law distribution with an universal 
index p ~ 2.2 — 2.3, we have determined the values of p from y-ray and X-ray spectra for 
a number of relativistic sources such as GRBs (prompt emissions and afterglows), blazars 
and pulsar wind nebulae. 

The maximum likelihood estimate of the width of the parent distribution for GRB prompt 
emission is found to be quite broad, dp = 0.51 ± 0.02; the probability that the distribution 
is consistent with a ^-function is extremely small, and therefore this result does not support 
that there is an universal p. 

We have considered the systematic errors in photon index due to spectra £t and time 
averaging of spectra and their contributions to the scatter in p distribution. We have shown 
that those contributions are very small for GRBs and can not explain the scatter in p 
distribution. 

For X-ray afterglows of GRBs, the p-distribution of the BeppoSAX sample can not rule 
out a possibility that the parent distribution is a 5-function distribution; however, a larger 
sample of Swift afterglows is inconsistent with a 5-function parent distribution. We point 
out that the smaller width of parent distribution for the BeppoSAX sample is due to its 
larger measurement error in j3. 

Analysis of 44 blazar spectra and 9 pulsar wind nebulae shows that the distributions of p 
for blazars and pulsar wind nebulae (PWNe) are also broad, not consistent with a 5-function 
distribution. 

Possible situations in which the “universality” of p could break are: (i) The shock is Mildly 
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relativistic (cf. Kirk et al. 2000); (ii) The magnetic field is oblique to the shock normal 
(Baring 2005); (iii) The nature and strength of the downstream magnetic turbulence are 
varying (Ostrowski & Bednarz 2002, Niemiec & Ostrowski 2004). A non-Fermi acceleration 
in a collisionless plasma shock was studied by Hededal et al. (2004), in which electrons 
are accelerated and decelerated instantaneously and locally, by the electric and magnetic 
fields of the current channels formed through the Weibel two-stream instability. It is not 
known whether an “universality” of p could hold for this mechanism. The “universality” 
of p might not happen in non-shock accelerations; for instance, in an alternative model for 
GRBs (Lyutikov & Blandford 2003), the energy is carried outward via magnetic field or 
Poynting flux. The particles accounting for the y-ray emissions are accelerated by magnetic 
field reconnection which may also produce a power-law spectra of accelerated particles with 
a variable p (however, this is still poorly understood). 
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12 3 4 


P 

Figure 1. Distributions of p for a sample of 394 GRB spectra with 100 < E^^ak < 200 keV. Solid line: using the relation 
p = —2(3 — 2. Dashed line: after narrowing the distribution by using the relation of either p = —2(3 — 2 or p = —2(3 — 1. 



100 < Ep^ak 

< 200 keV 

200 < Ep^ak 

< 300 keV 

Spectra samples 

Non-evolving 

Evolving 

Non-evolving 

Evolving 

< p > 

2.86 ± 0.06 

2.38 ± 0.03 

2.50 ± 0.07 

2.14 ± 0.03 

a(p} 

0.44 ± 0.04 

0.47 ± 0.03 

0.58 ± 0.06 

0.42 ± 0.03 


Table 1. Parameters of parent distribution of p for BATSE GRB spectra samples with Epeafe-^volution (AEp^ak > 15%E'pea/c) 
and without Epeafe-evolution (AEp^ak < 15%Epeafc)! where AEp^ak is the Ep^ak difference between any two adjacent-in-time 
spectra. All spectra are assumed in “fast cooling” regime. 
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Figure 2. The discrepancy between the fitted value and the ‘true’ value of /3, as a function of the higher end of the fitting 
energy range for the synchrotron spectra fitted by the ‘Band’ function. Two vertical lines mark the range of Emax!E-peak 
corresponding to the Ep^ak range of the sample in Fig. 1. The errors of the spectrum data are assumed to be proportional to 
square root of photon counts: <j{N{u)) oc N[v). 


1 - 1 - 1 - 1 - 1 -r 

(p)= 2.19±0.09 
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Figure 3. The distribution of p determined from 32 time-integrated GRB spectra. Solid line: p is inferred from the high-energy 
power-law index /3 by the relation p = —2(3 — 2. Dashed line: the narrowest distribution of p using the relation either p = —2(3 — 2 
or p = —2(3 — 1. /3 is taken from the ‘Band’-function fit by Band et al. (1993) to the time-integrated spectrum for each burst. 
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P P 

Figure 4. Left: The distributions of p for 11 GRBs observed by BeppoSAX (Amati et al. 2002); Right: The distributions of p 
for 13 X-ray flashes, X-ray-rich GRBs and GRBs observed by HETE-2 (Sakamoto et al. 2005). Solid lines: p is inferred from 
the higher-energy photon index 0 by the relation p = —2/3 — 2. Dashed lines: the narrowest distributions of p using the relation 
either p = —2/3 — 2 or p = —2/3 — 1. 




P P 

Figure 5. The distributions of p for GRB X-ray afterglows. Left: 14 afterglows are observed by BeppoSAX, taken from De 
Pasquale (2005). Right: 28 afterglows are observed by Swift, taken from O’Brien et al. (2006). Solid lines: p is inferred from the 
photon index /3 by the relation p = —2/3 — 2. Dashed lines: the narrowest distribution of p using the relation either p = —2/3 — 2 
or p = —2/3 — 1. 
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Figure 6. The distribution of p for 44 X-ray spectra of 33 blazars. Solid line: p is inferred from the photon index jS by the 
relation p = —2/3 — 2. Dashed line: the narrowest distribution of p using the relation either p = —2/3 — 2 or p = —2/3 — 1. /3 is 
taken from the catalog compiled by Donato et al. (2005). 
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